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SUMMARY 


This  report  covers  the  period  1  July  1971  -  30  June 
1972  which  is  characterized  by  system  developments 
and  software  debugging  efforts  and  establishment  of 
organizational  procedures  required  for  the  routine 
NORSAR  operation.  The  research  efforts  were  aimed 
at  system  improvements,  especially  to  increase  the 
event  detection  capability,  and  array  performance 
evaluation . 

Work  completed  and  in  progress  is  discussed  in 
Chapter  II.  The  first  five  sections  deal  with  event 
detection  problems,  while  in  the  next  four  sections 
observed  and  predicted  NORSAR  event  detection  and 
location  capabilities  are  discussed.  Two  studies 
on  event  classification  problems  are  also  included, 
and  the  last  four  sections  cover  structural  investi¬ 
gations  and  long  term  noise  variation  investigations. 

The  most  important  single  event  in  the  reporting 
period  was  the  seminar  on  Seismology  and  Seismic 
Arrays  which  took  place  in  Oslo,  22-25  Nov  1971. 


INTRODUCTION 

The  report  summarizes  the  NTNF/NORSAR  research  and  develop¬ 
ment  efforts  during  the  interval  1  July  1971  to  30  June 
1972.  In  the  first  part  of  the  period  most  attention  was 
given  to  development  work  like  software  debugging  and  modi¬ 
fication,  and  improving  the  operational  procedures  for  the 
array.  For  example,  a  computerized  log  of  the  operational 
quality  of  the  individual  SP  and  LP  seismometers  has  been 
created,  and  is  a  valuable  supplement  to  data  channel 
quality  coding  on  the  high  rate  tapes.  Recording  of 
NORSAR  LP  data  on  a  special  low  rate  tape  became  part  of  the 
real  time  operational  system  in  November  1971.  Similarly, 
on-line  data  transmission  from  SAAC  to  NDPC,  comprising 
processed  LASA  SP  and  ALPA  and  LASA  LP  data,  are  stored 
on  the  new  NORSAR  low  rate  tape.  A  complementary  event 
detector  is  under  implementation  in  the  Detection  Processor. 
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Several  modifications  of  the  Event  Processor  (EP)  have 
reduced  the  computer  time  per  processed  event  from  around 
9  min.  to  around  3.5  min.  v/ithout  appreciable  loss  in  the 
EP  performance.  As  a  result  of  the  above  work,  the 
NORSAR  operation  became  much  more  troublefree  in  1972, 
thus  providing  much  more  time  for  research. 


The  research  topics  investigated  or  in  progress  are  mainly 
aimed  at  system  improvements  and  evaluation  of  the  array's 
event  detection  and  location  performance.  Also,  some 
aspects  of  the  event  classification  problem  have  been 
considered.  Most  attention  has  been  given  to  improving 
NORSAR' s  event  detection  capability,  and  the  work  here 
comprised  theoretical  design  of  different  types  of  seismic 
event  detectors,  parameterization  of  detectors  for  best- 
in-average  performance,  event  detector  false  alarm  rate 
fluctuations  and  potential  SNR  gains  by  using  regionalized 
bandpass  filters.  The  so-called  incoherent  event  detec¬ 
tor  developed  by  NTNF/NORSAR  personnel  has  proved  useful 
for  real-time  multiarray  processing  problems,  as  demon¬ 
strated  by  analysis  of  three  events  recorded  by  WWSSN 
stations.  Procedures  for  simulating  array  event  detec¬ 
tion  and  location  capabilities  have  been  worked  out, 
including  an  off-line  version  of  the  on-line  Detection 
Processor.  A  comparison  between  observed  and  predicted 
NORSAR  event  detection  and  location  performance  gave 
essentially  identical  results.  Preliminary  results  on 
mj3:Ms  classification  capabilities  for  the  NORSAR  array 
have  been  obtained,  and  work  on  short  period  classifica¬ 
tion  criteria  has  started.  An  analysis  of  single  sensor 
amplitude  variations  across  NORSAR  favors  a  lognormal 
probability  density  distribution  function.  Vespagram 
analysis  of  core  precursor  phases  indicates  that  these 
waves  are  inadequately  explained  by  the  standard  models 
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for  the  earth's  core.  Mislocation  vectors  observed  for 
Fennoscandian  arrays  have  been  used  for  deducting  an 
upper  mantle  model  for  the  mentioned  region.  Also,  some 
results  on  long  term  noise  level  variations  have  been 
obtained . 

The  above  research  topics  and  relevant  results  are  described 
in  the  next  chapter.  In  general  only  the  main  results  are 
presented  here,  as  further  details  are  available  in  NORSAR 
Technical  Reports  or  from  papers  published  in  professional 
journals . 
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2 .  RESEARCH  EFFORTS 

The  research  activities  in  the  reporting  period  1  July 
1971  -  30  June  1972  have  been  focused  on  event  detection, 
localization  and  classification  problems.  The  NORSAR  event 
detection  performance  was  relatively  poor  in  the  fall  1971, 
mostly  due  to  high  noise  level  during  the  seasonal  storms 
in  the  North  Atlantic,  so  very  much  attention  was  given 
to  the  problem  of  improving  the  array's  detectability. 


The  mathematical  formulation  of  the  detection  algorithm 
in  use  is  as  follows: 

t 

STA(t)  -l  |  S  (i)  | 

i=t-IW+l 


where  t  =  STA  sampling  time,  S(i)  =  array  beam  amplitude 
and  IW  =  integration  window  length. 


LTA(t')  =  (1-2  n) -LTAft'-IWj+STAft'-IW) 


where  t'  =  LTA  sampling  time,  and  the  parameter  r)  =  5  out¬ 
side  detection  state  (no  signal  present) ,  and  n  =  4  in 
detection  state.  The  ratio  STA/LTA  is  calculated  every 
time  STA  is  updated.  When  this  ratio  exceeds  the  speci¬ 
fied  detection  threshold  a  certain  number  Q  of  consecu¬ 
tive  times,  a  detection  is  declared.  Currently  the  values 
of  IW,  STA  sampling  rate  are  15.5  and  1  respectively.  In 
view  of  general  detection  theory,  it  was  found  that  the 
optimal  detector  for  the  NORSAR  system  according  to  the 
Neyman-Pearson  criterion  is  a  square  detector  where  a 
whitening  filter  is  applied  on  the  beam  traces.  However, 
from  analysis  of  real  and  simulated  seismic  data  we  found 
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that  the  gain  was  very  modest  by  changing  from  a  linear 
to  a  square  detector,  and  the  same  applied  to  using  an 
extra  whitening  filter  in  the  NORSAR  system.  For  further 
details  on  this  problem  v/e  refer  to  a  thesis  by  Berteussen 
(1972)  . 


The_False_Alarm_Rate_of _the_NORSAR_Detector 

The  false  alarm  rate  at  NORSAR  does  vary  considerably  v/ith 
background  noise  level,  and  thus  reflects  fluctuations  in 
the  character  of  the  noise.  Otherwise,  if  the  noise  level 
increases  with  no  change  of  spectrum,  the  false  alarm  rate 
will  be  unchanged  as  long  as  the  detection  threshold  re¬ 
mains  fixed.  The  above  problem  was  investigated  by  R.T. 
Lacoss,  Lin.  Lab.,  M.I.T.,  during  a  visit  to  NORSAR  DPC 
in  April-May  1972.  An  idealized  conceptual  model  has  been 
formulated  to  help  understand  these  false  alarm  fluctua¬ 
tions  and  to  predict  certain  trade-offs  (Lacoss  1972) . 

For  example,  the  false  alarm  rate  under  best  noise  condi¬ 
tions  with  a  bandpass  filter  corner  frequency  of  0.9  Hz 
should  be  the  same  as  for  worst  noise  conditions  and  a 
filter  having  the  lower  corner  frequency  at  1.2  Hz. 


^_New_NORSAR_Event_Detector 

The  P-signals  recorded  by  NORSAR  are  only  partially  coherent 
across  the  array.  This  means  that  the  expected  gain  in 
signal-to-noise  ratio  (SNR)  which  is  proportional  to  the 
square  root  of  number  of  sensors  used  is  not  obtained 
during  array  beamforming  operations.  The  corresponding 
signal  energy  loss  increases  with  increasing  frequency, 
and  may  severely  degrade  the  array's  detectability  of 
very  short  period  P-waves.  This  problem  may  be  partly 
circumvented  by  replacing  or  supplementing  the  array  beam 
traces  with  the  average  of  subarray  beam  traces.  The  rela¬ 
tive  advantages  of  using  the  so-called  incoherent  beams 
are  modest  signal  losses,  better  estimates  of  the  noise 
variance  and  good  areal  coverage.  The  noise  suppression 


6 


is  small  as  compared  to  array  beamforming,  but  could  partly 
be  compensated  for  by  using  high  frequency  bandpass  filter¬ 
ing.  As  mentioned  previously,  an  experimental  event  detec¬ 
tor  based  on  incoherent  beams  was  implemented  in  the  on-line 
system  in  Dec  1971.  Preliminary  results  indicate  that  this 
new  supplementary  detector  when  completed,  would  improve 
the  array's  event  detection  performance  by  10-20  per  cent. 

A  theoretical  description  of  the  incoherent  detector  and 
its  operational  performance  has  been  prepared  by  Ringdal, 
Husebye  and  Dahle  (1972). 


In  addition  to  the  array  SNR  gain  resulting  from  beamforming 
processing,  additional  noise  suppression  is  obtainable  by 
applying  various  types  of  filters.  The  performance  of  a 
large  number  of  Wiener  prediction  error  filters  and  Butter- 
worth  3rd  order  bandpass  filters  have  been  checked  on 
several  types  of  frequency  domain  synthesized  noise  and 
P-wave  models.  In  the  latter  case,  the  signal  spectrum 
shape  is  fixed  while  the  peak  frequency  is  varied  in  fre¬ 
quency  steps  of  0.1  Hz. 

The  gain  in  SNR  due  to  prediction  error  (PE)  filtering 
is  not  very  impressive  except  when  the  signal  and  noise 
are  well  separated  in  the  frequency  domain.  The  reason 
for  this  is  simply  that  the  above  filter  is  essentially 
a  noise-whitening  filter.  A  single  bandpass  filter  gives 
in  general  a  somewhat  better  gain  in  SNR  than  a  PE  filter, 
but  a  stack  of  such  filters  is  required  for  optimum  per¬ 
formance.  For  example,  the  best-in-average  filter  has  a 
passband  of  ca.  1.4-3. 4  Hz,  but  might  be  5  dB  down  in  gain 
in  comparison  to  other  bandpass  filters  for  extremely  high 
and  low  frequency  P-signals,  like  those  from  Central  Ameri¬ 
ca  and  Central  Asia,  respectively. 

Using  the  NORSAR  observed  distribution  of  the  dominant 
P-signal  frequency  for  recorded  events  in  combination 
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with  relative  filter  gain  in  dB,  we  found  that  the  array's 
event  detectability  may  increase  with  roughly  10-20%  if 
"optimum"  bandpass  filtering  is  used.  The  above  work  was 
undertaken  by  Gj0ystdal  and  Husebye  (1972). 


Mul tiarraY_Processing_Experiment 

As  the  longitudinal  seismic  waves  are  not  dispersive,  the 
group  and  phase  velocities  are  identical  for  the  P  signals. 
This  means  that  under  certain  simple  assumptions,  i.e., 
the  case  in  which  the  amplitude  distribution  of  the  compo¬ 
nent  waves  in  the  wave  packet  is  Gaussian,  the  envelope  y 
of  the  signal  or  disturbance  is  expressible  as  (Coulson 
1961) : 

y  =  A(ir/a)2  exp  (-it2  (x-vt)  2/a) 

when  A  =  amplitude  scaling  factor,  a  =  STD  of  the  Gaussian 
distribution  and  x  is  distance.  The  half-width  of  the 

j. 

envelope  is  (a’/11)  and  thus  is  a  function  only  of  the  P- 
wave  spectral  bandwidth.  In  other  words,  we  would  expect 
the  signal  envelope  to  be  very  similar  from  one  station 
to  another,  even  if  the  signal  shape  may  vary  considerably. 

If  this  is  true  for  seismic  waves,  then  multiarray  or  global 
network  processing  would  be  feasible  using  an  event  detec¬ 
tor  similar  to  that  of  incoherent  beamforming  in  use  at 
NORSAR  as  discussed  in  a  previous  section. 

The  above  hypothesis  has  been  tested  on  WWSSN  station  records 
(two  earthquakes  and  one  explosion)  and  NORSAR  subarray  beams 
from  many  different  events.  Signal  envelope  similarity  has 
been  detected  through  cross-correlation  and  coherency  analy¬ 
sis,  and  a  typical  value  w as  around  0.75  units  between  WWSSN 
station  data.  Similar  results  were  obtained  by  analysis  of 
22  different  NORSAR  events  in  the  distance  range  3-145  deg. 
The  response  pattern  or  global  array  beampacking  was  also 
calculated.  In  short,  several  kinds  of  multiarray  pro¬ 
cessing  are  feasible  if  the  P-signals  are  replaced  by 
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their  envelopes.  The  results  outlined  above  we re  obtained 
by  Husebye,  Ringdal  and  Fyen  (1972) . 


One-ArraY_and_Two-ArraY_Location_Ca]Dabilities 

Event  location  capabilities  of  a  single  array  and  of  two 
arrays  in  terms  of  random  and  biased  observational  errors 
have  been  simulated  on  the  computer.  The  approach  v/as  to 
obtain  a  Gaussian  distribution  of  the  parameters  azimuth, 
slowness  (DT/DA)  and  arrival  time  differences  which  are 
used  in  one  and  two  array  epicenter  determinations.  Spe¬ 
cifying  their  mean  values  and  variances,  a  large  number 
of  paired  values  of  these  parameters  are  simulated  by 
using  a  random  number  generation  routine.  These  observa¬ 
tions  determine  a  distribution  of  event  locations  in 
geographic  space,  from  which  the  axes  of  the  95%  confi¬ 
dence  ellipse  were  calculated.  In  case  of  one  array,  the 
semi-axes  of  the  confidence  ellipse  are  tabulated  as  a 
function  of  epicentral  distance,  and  fictive  variances 
of  observed  azimuth  and  slowness.  The  joint  event  loca¬ 
tion  capability  for  two  arrays  is  easily  computed  but 
not  tabulated  due  to  the  regional  changes  in  the  shape 
and  orientation  of  the  confidence  ellipses.  The  above 
analysis  was  supplemented  with  epicenter  determinations 
using  real  NORSAR  and  LASA  data.  In  the  West  Indies  area 
the  average  location  error  was  found  to  be  310,  161  and 
139  km  for  NORSAR,  LASA,  and  NORSAR  and  LASA  combined. 

A  model  for  individual  subarray  time  corrections  for  the 
NORSAR  array  is  used  for  simulating  biased  correction  er¬ 
rors,  which  would  be  very  large  if  left  uncorrected.  The 
results  obtained  roughly  indicate  that  the  array  siting 
area  may  account  for  about  45  per  cent  of  the  observed 
slowness  anomalies,  but  no  improvement  was  obtained  in 
the  azimuth  anomalies.  The  slowness  bias  makes  it  dif¬ 
ficult  to  utilize  the  NORSAR  DT/DA  measurements  for 
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investigations  of  heterogeneities  in  the  mantle.  For 
further  details  on  this  problem  we  refer  to  a  forthcoming 
paper  by  Gj0ystdal  et  al  (1973)  . 


5YSiySti2Q_2^_£i22_§Y§12t_Detection_and_Location_Capabili  ty 
of _ the_NORSAR_array 

The  routinely  analyzed  NORSAR  data  from  the  first  G  months 
of  regular  operation,  May-October  1971,  has  been  evaluated. 

In  this  period  about  8  events  per  day  have  been  reported 
from  the  30°-90°  distance  interval.  This  gives  a  90%  cumu¬ 
lative  reporting  threshold  of  about  4.1  NORSAR  magnitudes, 
which  is  equivalent  to  aoout  4.3  NOAA  magnitudes.  The 
large  variations  in  the  seismic  background  noise  cause 
large  variations  in  the  distribution  of  reported  events 
per  time  unit.  Out  of  the  union  of  the  events  reported 
by  NORSAR  and  NOAA  in  May-August  1971,  25%  are  reported 
by  both,  while  30%  are  reported  only  by  NORSAR,  45%  only 
by  NOAA.  A  location  comparison  shows  that  50%  of  the 
NORSAR  locations  are  within  200  km  of  the  NOAA  location 
in  the  30°-90°  distance  range.  The  biased  errors  in 
time  delay  and  location  estimates  are  demonstrated  to  be 
quite  large,  which  indicates  that  significant  improvements 
in  the  performance  of  the  array  can  be  expected  when  a 
better  set  of  corrections  for  these  errors  is  implemented. 

The  above  investigation  which  is  based  on  NORSAR  data  from 
1971,  was  performed  by  Bungum  and  Berteussen  (1972)  .  How¬ 
ever,  during  January  1972  a  number  of  modifications  and 
improvements  were  implemented  in  the  Detection  and  Event 
Processor  software  systems  thus  necessitating  a  re-evaluation 
of  the  array's  operational  performance.  The  work  was  under¬ 
taken  by  Bungum  (1972a)  and  an  outline  of  the  main  results 
is  as  follows:  Data  published  in  the  NORSAR  seismic  bulle¬ 
tin  between  February  and  June  1972  has  been  studied  in 
order  to  find  estimates  of  the  detectability  and  location 
accuracy  at  NORSAR.  The  detectability  is  calculated  from 
empirical  frequency-magnitude  distributions,  and  the  90% 
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cumulative  detectability  for  the  teleseism.ic  zone  (30°<A<90°) 

has  been  estimated  at  m.  4.0,  while  values  for  different 

b 

regions  vary  from  3.7  to  4.3.  These  are  all  NORSAR  magni¬ 
tudes.  The  magnitude  bias  between  NOAA  and  NORSAR  has  been 
found  to  be  0.15+0.31  in  the  teleseismic  zone.  In  this 
zone,  the  median  location  difference  between  NOAA  and 
NORSAR  has  been  estimated  at  160  km,  with  values  for  dif¬ 
ferent  regions  ranging  from  130  to  340  km. 


E^§^i£t±2D_2^_tl}§_NQRSAR_Array  _Event_Det  ec  t  ability 

Obviously,  the  event  detection  capability  of  an  ordinary 
station  or  a  seismic  array  is  mainly  governed  by  the  noise 
level  at  the  site  under  the  assumption  that  instrument 
magnification  is  sufficiently  high.  In  NORSAR' s  case, 
factors  of  importance  are  besides  the  noise  level  on  the 
beam  traces,  the  threshold  to  be  exceeded  before  an  event 
detection  is  declared  and  the  relative  signal  losses  in¬ 
hibited  in  detection  processing  or  NORSAR  seismic  surveil¬ 
lance.  The  starting  point  was  to  estimate  the  cumulative 
distribution  of  the  linear  noise  power  parameter  STA  and 
its  relation  to  ground  motion  and  henceforth  event  magni¬ 
tude.  The  threshold  parameter  is  no  problem  as  it  is  a 
constant  while  an  estimate  of  the  relative  signal  loss 
was  obtained  by  comparing  Detection  and  Event  Processor 
STA  estimates  for  800  recorded  events.  A  comparison  be¬ 
tween  observed  (see  Bungum  1972  a)  and  predicted  events  for 
several  seismic  zones  are  shown  in  Table  1,  which  is  taken 
from  a  paper  by  Berteussen  and  Husebye  (1972) . 
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No. 

zonf: 
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ZONE  LI* 
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2 

At  1  antic 

180-260 

30-00 

38 

3.64 

4.76 

13 

3.69 

4  .23 

3 

N . Arcr iea 

260-340 

40-00 

114 

3.72 

4  .  qr 

! 

93 

3.66 

4.05 

4 

A  1  on  f  i  r*n  I‘i. 

340-15 

30-00 

131 

3.4" 

3.00 

17 

3.62 

3.95 

5 

Jaonr 

1  5-70 

50-00 

733 

3.^ 

4.07 

236 

3.61 

3.95 

G 

C.  Asia 

40-110 

30-00 

211 

3.21 

3.6° 

58 

3.45 

3.87 

7 

Iran  ( l ) 

110-160 

30-00 

26  2 

3.4  5 

3.80 

38 

3.51 

3.83 

8 

Iran  (2) 

110-130 

35-50 

!  133 

3.42 

3.78 

31 

3.49 

3.83 

TABLE  1 

Observed  and  predicted  detectability  levels  for  the  NORSAR  array. 
The  observational  data  (see  Bungum  1972  a)  covers  the  interval  Feb- 
June  1972. 


^YSDt_Classification_ApglYin2_the_mb-Kc._Criterion_to  NORSAR 

* 

A  study  has  been  made  by  Filson  and  Bungum  (1972)  of  the 
capability  of  the  Norwegian  Seismic  Array  to  discriminate 
between  earthquakes  and  underground  explosions  occurring 
in  Central  Asia  and  western  Russia.  The  ratio  of  surface 
to  body  wave  magnitudes  (Ms:m  )  has  been  used  exclusively, 
and  measurement  of  low  amplitude  surface  waves.  Beamforming 
and  matched  filtering  were  the  signal  enhancement  techniques 
applied.  Of  34  events  in  central  Asia  studied,  10  were  iden¬ 
tified  as  explosions,  22  earthquakes  and  2  were  unidenti¬ 
fiable  because  of  high,  long  period  background  noise. 

Five  events  in  aseismic  western  Russia,  all  presumed  to 

be  explosions,  showed  wide  variation  in  M  :m,  two  of 

s  b 

these  measurements  being  close  to  the  earthquake  population 
of  central  Asia. 


*  J.  Filson,  Lincoln  Lab.,  M.I.T.,  paid  a  visit  to  NORSAR 
Data  Processing  Center  Oct-Nov  1971. 
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The  Ms  rm^  measurements  of  central  Asian  explosions  made 

in  this  and  another  study  are  compared  v/ith  measurements 

of  Nevada  explosions.  For  a  given  m^  the  latter  consistently 

have  hicTher  M  values.  The  possible  causes  of  this  are  dis- 
s 

cussed  and,  depending  on  the  cause,  the  array  seems  capable 
of  detecting  surface  v/aves  from  an  explosion  of  4-16  kilo- 
tons  yield  at  40°  distance  in  central  Asia.  Strong  seis¬ 
mic  noise  variations  in  Norway  make  this  long  period  capa¬ 
bility  time  dependent. 


§i}2^t_Period_Event  _Uiscrim.inants 

The  design  of  short  period  seismic  event  discriminants  is 
principally  aimed  at  manifesting  typical  source  charac¬ 
teristics  of  underground  nuclear  explosions  and  earthquakes. 
A  special  problem  here  and  typical  for  an  array  like  NORSAR 
where  the  P-waves  are  only  partially  coherent,  is  the  signal 
energy  loss  during  ordinary  data  processing.  For  example, 
high  frequency  signal  energy  is  strongly  suppressed  during 
conventional  array  beamforming,  but  preserved  on  the  sub¬ 
array  beam  level.  Henceforth,  our  initial  event  discrimi¬ 
nant  efforts  are  focused  on  processing  which  minimized 
the  loss  of  the  important  high  frequency  components  in  P- 
waves .  Using  the  Fast  Fourier  Transform  as  a  basic  tool, 
power  spectra  of  short  period  P-wave  signals  from  16  well- 
recorded  earthquakes  and  underground  explosions  in  Central 
Asia  and  in  western  North  America  were  computed.  The  ini¬ 
tial  data  preparation  included  a  proper  line-up  of  incoming 
signals  and  tapering  the  ends  of  the  6.4  sec  transform 
windows  used.  The  spectra  were  smoothed  and  the  effect 
of  the  noise  carefully  removed.  The  computed  spectra  are 
the  mean  power  of  all  132  sensors  in  the  NORSAR  array 
(spectraform)  and  the  power  on  the  array  beam.  The  ratio 
of  the  above  two  spectra  (beam  loss)  increases  sharply  at 
frequencies  above  1.5  Hz  being  in  average  12  dB  at  3  Hz 
Earthquakes  and  explosions  were  observed  to  have,  in 
average,  different  shapes  of  spectrum.  The  observed 
spectraform  spectra  peaked  between  1  and  2  Hz.  The 
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signal-to-noise  ratio  for  the  events  in  Central  Asia 
peaked  around  2  Hz,  but  was  more  flat  for  events  in 
North  America  due  to  shift  of  energy  towards  lower  fre¬ 
quencies.  At  frequencies  above  the  spectral  peak,  spectra- 
forming  gave  the  most  stable  estimates,  though  a  combined 
processing  method,  the  average  spectrum  of  subarray  beams, 
also  had  high  stability.  The  above  work  has  been  per¬ 
formed  by  Noponen*,  Husebye  and  Rieber-Mohn  (19  72)  . 


Ti™§_32l3_Ami2litude_Anomal  ies_on_the_NORSAR_SubarraY_Level 

An  investigation  is  undertaken  to  map  time  delays  and 
amplitude  anomalies  on  the  subarray  level,  and  this  means 
that  all  132  NORSAR  SP  sensors'  performances  in  this 
respect  are  jointly  analyzed  for  the  first  time.  The 
main  results  obtained  so  far  are,  according  to  A.  Dahle, 
as  follows.  The  single  sensor  amplitude  varies  tremen¬ 
dously  across  the  array,  a  factor  of  2  and  10  on  the 
subarray  and  array  level  is  not  unusual.  However,  the 
interesting  feature  is  that  Kolmogorov-Smirnov  hypo¬ 
thesis  tests  favor  a  lognormal  probability  distribution 
for  the  observed  amplitude,  and  this  is  very  clear  for 
high-frequency  signals  (see  also  Ringdal  et  al  1972) . 

Time  anomalies  across  the  subarrays  defined  as  the  dif¬ 
ference  between  observed  and  predicted  from  Herrin's 
tables  are  in  general  small,  i.e.,  the  standard  devia¬ 
tions  were  less  than  around  0.1  sec.  The  signal  amplitude 
loss  due  to  subarray  beamforming  is  in  average  less  than 
1  dB  and  varies  linearly  as  a  function  of  the  mean  cross¬ 
correlation  value. 


*  I.  Noponen,  Helsinki  University,  visited  NORSAR  Data 
Processing  Center  in  the  intervals  16  Aug  -  17  Dec  1971 
and  20  May  -  30  June  1972. 
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Vcsgagram_AnalYsis_of _Core  Waves 

Some  problems  in  the  interpretation  of  core  phases  and 
their  precursors  have  been  studied  by  the  analysis  of 
PKP  phases,  recorded  at  the  NORSAR  array.  A  limited 
number  of  SKP  phases  has  been  used  as  an  aid  in  the 
analysis.  Processing  techniques,  suitable  for  the  reso¬ 
lution  of  core  phases  in  slowness  and  time  are  discussed, 
in  particular  the  usefulness  and  limitations  of  the  Vespa 
process  in  this  connection.  In  the  various  techniques 
the  effect  of  the  structure  under  the  array  should  be  taken 
into  account.  The  methods  used  reveal  the  refractions  and 
reflections  at  the  inner  core  boundary  and  also  to  some 
extent  the  more  complicated  character  of  the  PKP  precur¬ 
sors.  Several  representative  examples  of  interpretation 
of  the  precursors  are  viewed  in  light  of  the  observational 
results.  The  analysis  is  limited  to  A>136°.  In  this  range 
the  results  indicate  that  diffraction  effects  may  play  a 
more  important  role  in  the  generation  of  the  precursors 
than  generally  assumed.  For  more  details,  we  refer  to 
a  paper  by  Doornbos*  and  Husebye  (1972) . 


yBE§£_^52i'le_Discontinuities_in_Scandi.nav:i.a 

During  a  research  visit  to  NORSAR  DPC,  Dr.  I.  Noponen 
investigated  slowness  anomalies  observed  at  array  stations 
in  Scandinavia,  namely,  NORSAR,  Hagfors  and  the  Helsinki 
tripartite  station.  His  main  results  are  as  follows: 
(Noponen  1972)  :  Systematic  teleseismic  P-wave  slowness 
anomalies  up  to  1  sec/deg  are  observed  at  three  array 
stations  in  Fennoscandia,  NORSAR  in  Norway,  Hagfors  Ob¬ 
servatory  in  central  Sweden  and  the  telerecording  station 
at  Helsinki  in  Finland.  The  slowness  anomalies  are  smallest 


*  Dr.  D.  Doornbos,  Utrecht  University,  visited  NORSAR 
Data  Processing  Center  in  the  summers  of  1971  and 
1972. 
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at  Helsinki  and  largest  at  Hagfors .  Resulting  event  mis- 
locations  would  exceed  ten  degrees  at  certain  regions  if 
no  corrections  were  used.  It  is  shown  that  an  integration 
over  the  horizontal  velocity  gradient  along  the  ray  path 
near  the  receiver  gives  the  slowness  anomaly  caused  by 
that  gradient.  A  given  horizontal  gradient  produces  an 
order  of  magnitude  smaller  mislocation  when  situated 
near  the  source  than  when  situated  near  the  receiver.  It 
is  concluded  that  a  major  part  of  the  observed  slowness 
anomalies  can  be  explained  by  lateral  changes  in  structures 
in  the  upper  mantle  and  the  crust  under  the  arrays.  In 
particular,  an  increase  of  the  average  vertical  velocity 
in  the  upper  mantle  under  the  Scandinavian  peninsula  in 
the  direction  from  oceanic  towards  continental  areas  is 
required  by  the  slowness  anomalies  and  by  the  travel  time 
residuals  of  the  conventional  stations  in  the  same  region. 


In  a  recent  paper  Bungum  (1972  b)  discusses  different  types 
and  principles  for  microseismic  noise  investigations.  The 
emphasis  is  on  arrays  like  LASA  and  NORSAR  which  may  be 
considered  very  efficient  tools  for  analyzing  the  micro- 
seismic  noise  fields  in  the  f requency-wavenumber  domain. 
Another  interesting  topic  discussed  in  the  above  paper  is 
the  long  term  variation  of  the  NORSAR  array  beam  noise 
level,  which  has  a  direct  bearing  on  the  array's  event 
detectability  as  discussed  in  a  previous  section. 

Spectral  analysis  of  short  period  noise  level  variations 
in  the  passband  1.2-3. 2  Hz  gives  a  significant  peak  at 
1  c/day,  and  this  reflects  the  diurnal  noise  variation. 
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3  .  MISCELLANEOUS 

To  mark  the  formal  opening  of  the  NORSAR  array,  a  seminar 
on  Seismology  and  Seismic  Arrays  was  arranged  in  Oslo, 
22-25  November  1971.  A  total  of  69  scientists  from  12 
countries  attended  the  meeting  and  altogether  28  talks 
were  given.  Twenty- two  of  the  papers  presented  at  the 
seminar  have  been  edited  by  Husebye  and  Bungum  (1972)  in 
a  special  booklet  titled  Proceedings  from  the  Seminar 
on  Seismology  and  Seismic  Arrays,  Oslo,  21-25  November 
1971 . 
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